Class formation in a social network with asset exchange. 
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In this paper the study of two kinds of economic exchange, additive and multiplicative, in a system 
of N agents has been made. The work is divided in two parts, in the first one, the agents are free 
to interact with each other. The system evolves to a Boltzmann-Gibbs distribution with additive 
exchange a condenses with a multiplicative one. If bankruptcy is introduced, both types of exchange 
lead to condensation. Condensation times have been studied. In the second part, the agents are 
placed in a social network. We analyze the behavior of wealth distributions in time, and the 
formation of economic classes was observed for certain values of network connectivity. 
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I. INTRODUCTION 

We have had in the past years a large amount of litera- 
ture dealing with the study of the distribution of wealth 
in agent based models with various kinds of interaction 
rules ^ 12 E]- Several distributions such as Boltzmann 
- Gibbs, Gamma, or Pareto can be obtained according 
to the different conditions of the models [U |5J [B] . It is 
well known that real data analysis from several countries 
yield a Boltzmann-Gibbs distribution for that sector of 
the population with lowest wealth, who are the majority, 
and a Pareto distribution for the minority of the popu- 
lation with the highest wealths. This particular behav- 
ior has been reproduced, to some extent, using different 
kinds of assumptions ^ . 

In this article we describe a model where this behavior is 
obtained and the appearance of social classes is observed. 
These two phenomena are related as we shall see in the 
following. We analyze the conditions under which the 
different types of distributions are obtained and also the 
conditions necessary for the social classes to appear. In 
the next section of the paper we describe the interacting 
model for the several cases that we study. We focus on 
the behavior of entropy, poverty (defined as the minimum 
state of wealth, where the amount of money is less than 
minimum allowed exchange) and wealth distributions in 
sections III and IV. 

In section III the dynamics of the model with and without 
bankruptcy is investigated and in section IV we study 
the effect of a network in the system. Finally in the last 
section we give our conclusions. 
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II. INTERACTING MODEL 

A. Interacting model with no network 

The kind of interactions that we have considered in the 
present work are those in which the total amount of 
wealth before and after the interaction is conserved, also 
known as elastic collisions. That is, if are the labels 
of the two agents involved and their wealths are (wi,Wj) 
respectively, then we can write 

Wi{t + M) ^ w^{t) + /:^w, (1) 
Wj{t + At) = Wj{t) - Aw. (2) 

Since wealth is preserved in the interaction, we have that 
Wi{t + At) + Wj{t + At) = w^{t) + Wj{t). It is worth 
to mention that, according to which amount we choose 
for Aw, we can have two distinct processes: additive or 
multiplicative exchange. 

In the additive exchange we have Aw = c, this means 
that the exchange money is fixed in the whole process, 
which produces an interaction with time symmetry, that 
is, if one interaction is produced, for instance i wins and 
gets Aw from j, in the next time step it is possible to 
have the inverse one, obtaining the original situation. Of 
course, since we are dealing with a system with a large 
number of agents, the probability of that to occur is very 
low (l/4Af^). In the multiplicative case since Aw is given 
by Aw — round{i' ■ min(w, . ifjj))|25j. < v < 1, time 
symmetry is broken and the system never returns to its 
original set of values. 

B. Interacting model within a network 

We introduce a network in the system in which every 
agent i has Ui links to other agents, such that 1 < < 



2 



o-maxi where amax is the highest possible number of Hnks 
allowed in the dynamics. Since this number is not a con- 
stant, according to the way the links were introduced 
(randomly) , the average number of links (a) in the net- 
work is given by 



-E 



(3) 



With this in mind the dynamics in the network is of 
course different from case |II A| since agent i will only 
interact with her links through equations ([l]) and (|2|. 
The simulation begins with a population of N agents. 
Each one of them have the same initial amount of wealth 
(w). In a Monte Carlo step (MCS or time step), two 
agents (i,j) are chosen randomly and interact through a 
bet. In the bet one of the agents is chosen randomly to 
be the winner (of course the remaining agent is the loser) , 
and an amount of wealth Aw is transferred from the loser 
to the winner. If the wealth of the loser is less than Aw, 
the transaction is not done. It is important remark that 
agents are not allowed to have negative wealth. 
In all simulations we have used N = 500, (w) = 100 and 
lO'^ realizations during 4 x 10^ MCS. In the additive ex- 
change the value of c is 20, in the multiplicative exchange 
v = 0.2. Within a network the average number of links 
used were 1.468, 1.952, 2.428 and 10.06 (these values cor- 
respond to values of amax ~ 2,3,4 and 20, respectively). 



III. DYNAMICS WITHOUT NETWORK 



A. Exchange without bankruptcy 
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in the additive case the system evolves to a stationary 
state and the distribution is of the Boltzmann-Gibbs type. 
In the multiplicative case the system undergoes differ- 
ent phases including the Gamma distribution (Fig[T|(c)) 
and the power law (Fig 0(d)). Also, as al- 

ready noted in previous works |12j . the system ends in 
a condensed state, where N — 1 agents belong to the 
poverty state, and a single agent is in the maximum state 
of wealth. 

Wc have also studied the Shannon entropy, defined by 



5 = -^PfelnPfe, 



(4) 



where Pk = Uk/N , that is, the probability for the agent 
to be in the k state, with wealth w^. In Figure |2ja) and 
|2jb) we display the behavior of the Shannon entropy as 
well of the poverty for the additive and multiplicative ex- 
changes. In the multiplicative case the entropy reaches 
a maximum before decreasing monotonically, in conse- 
quence, the number of states available to the system is 
minimum at the condensed state, Sc |26j . In our case, 
with the parameters given for the system, this value is 
Sc = 0.01442 
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Figure 2: Time evolution of entropy (a) and poverty (b) for 
both cases of exchange (additive and multiplicative). 



B. Exchange with bankruptcy 
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Figure 1: Time evolution of wealth distribution for both cases 
of exchange (additive and multiplicative). 

We studied the wealth distribution evolution with respect 
to time for both the additive and multiplicative exchange 
cases. In Figure [T] we display the results, which are not 
normalized, and find out what it is already reported [TJ[H]: 



We say that an agent is bankrupt when it ends up with 
no money to make more bets; in that case, it is not al- 
lowed to participate in the dynamics. As expected, the 
condensation point is reached faster than in the previous 
cases when bankruptcy was not included. In those cases, 
the agents who have reached the absolute poverty were 
still allowed to participate in the game, and in principle 
they could become rich, even though this situation was 
very unlikely. 
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Figure 3: Time evolution of entropy (a) and poverty (b) for 
both cases of exchange (additive and multiplicative), using 
the bankruptcy concept. 
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Figure 5: (a) Time of condensation vs Aw (additive), (b) time 
of condensation vs i/ (multiplicative), using the bankruptcy 
concept. 



Figure [3] shows the results with bankruptcy. It is impor- 
tant to note that in both cases, additive and multiphca- 
tive, the system reaches the condensed phase. Poverty 
reaches its maximum value, — 1 or 499, and entropy 
ends up at the Sc value previously given. Since in both 
cases we obtain condensation, we consider the time it 
takes to reach this point when we vary parameters such 
as N, the number of agents, and Aw, the fraction of the 
wealth exchange. In Figure[4]we show this variation with 
respect to N and found a different behavior depending of 
the two cases. In the additive case the dependence of the 
condensation time with respect to N goes like t'^'^'^ oc 
whereas in the multiplicative case we get a linear depen- 
dence N f™"' cx N. The dependence of the condensation 
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Figure 4: Time of condensation vs number of agents for 
both cases of exchange (additive and multiplicative), using 
the bankruptcy concept. 

time with Aw is depicted in FigurcjSj In the additive case 



(fig[5ja)) there is a discrete behavior due to the fact that 
the condensation time has a well defined value for the dif- 
ferent ranges of Aw. This can be explained if we consider 
the condensation time as being proportional to the mini- 
mum of the exchange [27J the two agents can make before 
one of them reaches absolute poverty tc oc minex- For 
instance if we considered (w) — 100 and Aw = 50, two 
agents (?,j) will have a minimum of exchange mirier = 2, 
that means that there are 2 exchanges at least, before 
one agent reaches poverty. If we have (w) = 100 and 
Aw = 51, the minimum of exchange is now miriex = 1, 
and only one bet can be made. Therefore the time of 
condensation depends of the minimum of exchange. 
In a multiplicative exchange (Figure [5f^b)) the behavior 
of the condensation time is of the form i™"' oc 1 /v'^ and 
do not shows the discrete situation of the additive case 
since the bet is not constant. One would expect that if 
we increase the amount of the wealth exchange in each 
transaction, the time it takes to condensate gets lower, 
and that is what the above relation is telling us. 



IV. DYNAMICS OF THE EXCHANGE IN A 
SOCIAL NETWORK 

In this section we consider a system composed of N 
agents interacting with the same rules as before, but this 
time they cannot interact freely but instead they have 
relations that form a social network, that we are about 
to discuss. One would expect that the dynamics shows 
some differences with respect to the urn case. 
We build the network in the following way. A number 
of links is chosen randomly by every agent in the popu- 
lation, this number ranges from one up to a maximum 
value of links amax- Now the agents can do bets only 
with other agents to whom they are linked, generating 
the social network. 
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Since the number of links varies among the agents, we 
will consider an average value of links, (a) (see section 
for a given network. We consider again the additive 
and multiplicative cases, without taking bankruptcy into 
account, because as we have seen, its main effect is to 
favor the system to go into the condensed state. 



A. Additive exchange 
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Figure 6: Time evolution of wealth distribution in the additive 
exchange case for different values of (a). 

In Figure [6] we display the temporal evolution of the 
wealth distribution for the additive case and for differ- 
ent values of (a). Again the distributions are not nor- 
malized since this makes no difference in the discussion. 
We always get a stationary distribution of the wealth 
no matter the value of (a). The entropy, as well as the 
poverty, reaches a maximum value (Figure [7]). There IS 
also an interesting phenomenon which is the "isle" for- 
mation, which consist of groups of linked agents that are 
isolated from the rest of the network. The exchange be- 
comes isolated from the rest of the system, and the effect 
of this situation is to produce economic classes as the 
system evolves. (See inset of Figure [?]). This effect is 
more easily obtained when (a) is lower, as one would ex- 
pect. Therefore, the isolated agents that form the isle are 
decoupled from the system's dynamics, and they never 
reach absolute poverty. The more isles we get the less 
poverty we obtain. 

Figure [8] shows wealth distributions for additive exchange 
for different values of (a). For the lower value of (a) — 
1.468, we noticed two social classes, one where the wealth 
of the agents is below 2 (w) and with a Boltzmann - Gibbs 
distribution (continuous line Figure ^a)) of the form 



P{w)^Ce-'^, (5) 
where k is equal to 2. The other class, with wealth above 
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Figure 7: Time evolution of entropy (a) and poverty (b) in 
the additive exchange case for different values of (a). 
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Figure 8: Wealth distribution in the additive exchange case, 
in time step 399000 MCS for different values of (a). 



2 {w) behaves like a Pareto distribution (dashed line Fig- 
ure [SF a)) of the form 



P{w) 



C 



(6) 



with a — 2.5. These two distributions are separated 
at a value of wealth w — 2 (w) , and it is important to 
mention that these classes coexist. Other authors have 
found a transition from one distribution to another as 
the system evolves |7], but not at the same time. These 
formation of two classes arises from the fact that at a 
value of (a) — 1.468, most of the isles are comprised of 
two agents. Therefore, the poor sector of the system is 
formed mostly by isles of 2 agents, while the rich sector 
have to the possibility to be formed by isles of 3 or more 
agents, which allows an agent to have a wealth higher 
than 2 (w) . 
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When the isles formed by two agents start to disappear 
(figures [8|D-|8ji) , also the two classes do vanish. Further- 
more when the value of (a) changes, also the value of k 
(eq ([5])) does, as shown in table |l] 



(a) 


k 


1.952 


1.15 


2.428 


1.1 


10.06 


1.05 



Table I: Values of k for different average links in the additive 
exchange case 



Note that the systems reaches equilibrium when the 
"temperature" or k (w) is higher than in the non-network 
case, which has fc = 1. As (a) increases, the value of k 
tends to 1 , which makes sense since a very large value of 
(a) is equivalent to having no network present. 
It is important to remark that the formation of the isles 
in the network allows the exchange transactions to occur 
in a local way and that prevents a population of agents to 
become poor. The agents are in a Boltzmann - Gibbs dis- 
tribution, since in the isles the maximum possible wealth 
is 2 (w) . Therefore isle formation prevent agents from 
becoming poor, but also rich. 



the four parameters of (a) (1.468, 2.428 and 10.16). It 
is also important to note that these distributions are not 
stationary when the entropy is at its maximum value, 
but they are at later times at values shown in Figure |10[ 
In Figure 10 one can see that entropy goes through 



a maximum before decreasing and reaching a constant 
value that is different for each value of (a). The system 
with larger parameter (a) = 10.06 is the one with more 
possibilities to condense, as in the non-network case. As 



one can see from Figure 10 ^a) and 10 ^b), its entropy is 



the lowest and its poverty the highest 
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B. Multiplicative exchange 



Figure 10: Time evolution of entropy (a) and poverty (b) in 
the multiplicative exchange case for different values of (a). 
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Figure 9: Time evolution of wealth distribution in the multi- 
plicative exchange case for different values of (a). 

In Figure |9] we show the time evolution of the wealth dis- 
tribution in a multiplicative exchange process for differ- 
ent values of (a). One can observe that in the same way 
as in the additive case, the distributions are very similar 
in the first time steps, however, they soon differentiate 
from each other. For instance, for a time of 10000 MCS, 
the wealth distributions are different at least for three of 
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Figure 11: Wealth distribution in the multiplicative exchange 
case, in time step 399000 MCS for different values of (a). Note 
that {w) = 100. 

Figure [TT] shows the wealth distributions in a multiplica- 
tive exchange at the same time step {t = 399000 MCS) 
for different values of (a) . A very interesting phenomenon 
can be observed in Figure 1 1 'a) , the appearance of three 



well defined classes with wealth of zero (not visible in the 
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figure), 1.8 (w) and 2.8 (w) respectively (isolated points 
in the figure). These classes of wealth will eventually be- 
come 0, 2 {w) and 3 {w) when the system reaches the sta- 
tionary state at later times. The isolated point of wealth 
1.8 (w) corresponds to the presence of isles of two agents 
in the system, because, if two agents interact only with 
each other, the maximum value of wealth that a agent 
will have is 2 (w) . Similarly, the point with wealth equal 
to 2.8 {w) corresponds to the presence of closed groups of 
three agents, where an agent will have the possibility to 
get at most 3 (w) of money. In addition to these classes, 
we observe a background of agents with variable wealths. 
These background arises from agents that do not belong 
to closed groups. In this case there is a larger fiux of 
money among many agents in the network. 
It is clear that classes appear due to the formation of 
isles. A local condensation occurs in each isle, and this 
gives us a hint about the way the system behaves when 
the value (a) decreases: the entropy of the system should 
reach a minimum which is stationary and in this point 
half of the agents are poor and the other half has a wealth 
of 2 {w). As soon as (a) starts to increase, the separation 
between the classes begins to disappear, until we are left 
with a class of agents with wealth 0, and another class 
on non-zero wealth agents, as seen in Figure [Tl^d) . A 
Similar bimodal behavior has been obtained in Ref [S] . 



V. SUMMARY AND CONCLUSIONS 

In this work we studied an asset exchange model with 
additive or multiplicative exchanges. We consider the 
cases when bankruptcy is introduced and also the inter- 
action of agents through a social network. As already 
known, the additive, non-network case leads to a Boltz- 
mann - Gibbs distribution, while in the multiplicative, 
non-network case, the system evolves to a condensed 
state where a single agent gets all the money. When we 
allow for bankruptcy, condensation is obtained for both 
types of exchange, again in the non-network case. We 
studied these condensation times, finding that the in the 
multiplicative exchange the system condenses faster than 
in the additive case. The time of condensation is propor- 
tional to the number of agents. In the additive case the 
proportionality goes as i"'^'' oc ■ minex where, mitiex 
is the minimum exchange two agents can have before one 
of them reaches poverty. For the multiplicative case the 
condensation time goes as t"™' cx N/v'^ . 
From the analysis presented for the wealth distributions 
we have obtained information about their behavior for 
different values of time and average value of links, (a). 
We analyzed two different processes, the multiplicative 
and the additive cases, with and without social network 
among the interacting agents. When we introduced the 
social network we obtain stationary distributions in the 
additive case with a clear formation of coexisting classes, 
the poor sector of the system following a Boltzmann- 
Gibbs distribution until a value of wealth 2 {w) and a 



Pareto distribution from there on. 

In the multiplicative case, there is also formation of three 
well defined classes that arise from the formation of isles 
in the network, and also from groups of 3 or more agents 
isolated from the rest of the network. These isles and 
groups begin to disappear when the parameter (a) is in- 
creased, since the connectivity of the network is higher. 
In this later case we expect a global condensation to oc- 
cur a very high times. 

Isle formation in general allows the presence of social 
classes, since they prevent the money to fiow in the en- 
tire network. This has as a consequence that in some 
parts of the system the wealth can not be larger than 
a maximum amount fixed by the number of members in 
the isle, effectively preventing global condensation. The 
number of isles is higher as (a) diminishes. In a real 
society, the flow of money is affected by a myriad pro- 
cesses that cannot be accounted for in a model. All these 
processes together lead to the formation of classes in a 
real economy, however, it is important to model some 
of the phenomena, since this will give us some insights 
into the complex whole that make up the total behavior 
of economic activity. The model developed in this work 
has the virtue of producing an interesting phenomenon 
like the formation of economic classes, despite the fact 
that it is so simple and does not take into account the 
whole economic complexity. We have also observed a 
correlation between isles and the formation of classes, so 
we would expect that isle economies would be present in 
real economies; in other words, sectors of the population 
whose interaction with the economy is negligible. This 
point could be explored in a future work. 

Acknowledgments 

This work has been partially supported by CONACyT- 
Mexico under grant 45782 and PROMEP-Mexico under 
grant NPTC-256. 



[1] V. Yakovcnko (2007), arXiv:0709.3662v3 [physics.soc- 
ph]. 

[2] A. Gupta (2006), arXiv:physics/0604161v2 [physics.soc- 
ph]. 

[3] M. Rodn'gucz-Achach and R. Huerta-Quintanilla, Phys- 

ica A 361, 309 (2005). 
[4] A. Chattcrjco and B. Chakrabarti, Eur. Phys. J. B 60, 

135 (2007). 
[5] A. Gupta, Physica A 359, 634 (2006). 
[6] M. Patriarca, A. Chakraborti, and K. Kaski, Phys. Rev. 

E 70, 016104 (2004). 
[7] J. Gonzalcz-Estcvcz, M. Cosonza, R. Lopez-Ruiz, and 

J. Sanchez, Physica A 387, 4637 (2008). 
[8] A. Dragulescu and V. Yakovcnko (2003), arXivxond- 

mat/0211175vl [cond-mat.stat-mech] . 
[9] Y. Brenner, J. Reijnders, and et al, The theory of incom- 
ing and wealth distribution (Martin's Press, 1988). 
[10] Chatterjee, Arnab, Yarlagadda, Sudhakar, Chakraborti, 

and B. K, Econophysics of Wealth Distributions 

(Springer Milan, 2005). 
[11] M. Patriarca, A. Chakraborti, and K. Kaski, Physica A 

340, 334 (2004). 
[12] C. Mourkazel, S. Gongalves, J. Iglesias, M. Rodriguez- 

Achach, and R. Huerta-Quintanilla, Eur. Phys. J. Special 

Topics 143, 75 (2007). 
[13] S. Ispolatov, P. Krapivsky, and S. Redner, Eur. Phys. J. 

B 2, 267 (1998). 
[14] A. Drogulescu and V. Yakovenko, Eur. Phys. J. B 17, 

723 (2000). 

[15] P. Mohanty, Phys. Rev. E 74, 011117 (2006). 

[16] D. Garlaschelh and M. I. Loffredo, J. Phys. A 41, 224018 

(2008). 

[17] F. Rcif, Fundamentals of statistical and thermal physics 

(McGraw Hill, 1965). 
[18] L. Brillouin, Science and information theory (Academic 

Press, 1962). 

[19] D. Ycvick, A first course in computational physics and 
object-oriented programming with C++ (Cambridge Uni- 
versity Press, 2005). 

[20] C. Gardiner, Handbook of stochastic methods for physics, 
chemistry and the natural sciences (Springer Verlag, 
1997). 

[21] N. V. Kampcn, Stochastic processes in physics and chem- 
istry (Elsevier, 2007). 

[22] H. Gould, J. Tobochnik, and W. Christian, An intro- 
duction to computer simulation methods applications to 
physical systems (Addison Wesley, 2006). 

[23] J. Gentle, Random Number Generation and Monte Carlo 
Methods (Springer Verlag, 1998). 

[24] D. Knoke and S. Yang, Social Network Analysis (SAGE 
Publications, 2008). 

[25] We used round{x) as the integer closest to x. 

[26] The entropy of condensation Sc is calculated as: 

/AT-IAT-l 1,1\ 
[27] The minimum of exchange is given by 



